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We uncover a new type of unitary operation for quantum mechanics on the 
half-line which yields a transformation to "Hyperbolic phase space" {r],Pr,)- 
We show that this new unitary change of basis from the position x on the half 
line to the Hyperbolic momentum p^, transforms the wavef unction via a Mellin 
transform on to the critial line s = 1/2 — ip^. We utilise this new transform 
to find quantum wavefunctions whose Hyperbolic momentum representation 
approximate a class of higher transcendental functions, and in particular, ap- 
proximate the Riemann Zeta function. We finally give possible physical reali- 
sations to perform an indirect measurement of the Hyperbolic momentum of 
a quantum system on the half-line. 

1 Introduction 

Higher transcendental functions find many uses in mathematics, engineering, physics and many 
other sciences. Their efficient numerical evaluation is a highly challenging task and typically 
one makes use of pre-compiled library routines (i). In this work we present a curious method 
which can design a ID quantum system in such a way that the resulting wavefunction is propor- 
tional to a given transcendental function taken from a certain class of functions. The resulting 
quantum system displays enormously rich behaviour. By studying the corresponding quan- 
tum system one might discover new insights into the properties these transcendental functions. 
More interestingly we discover how to execute a new quantum unitary transform, the quantum 
Mellin transform, (similar to the quantum Fourier transform), which takes the Mellin trans- 
form of a wavefunction. The quantum Fourier transform has been prominent in many quantum 
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algorithms and one might suspect that the quantum MeUin transform might also be useful in 
quantum computation. 

The class of transcendental functions we consider are the three parameter family of higher 
transcendental functions sometimes known as the Lerch transcendents $(2;, s, u), which contain 
the celebrated Riemaim Zeta function ( (2) section 1.11 and (5) section 64:12). The Lerch 
transcendent is given as the analytic continuation of the series 

00 ^ 

^z,s,u)^y^ ^ ,kl<l, 1^7^ -1,-2,-3,- •• (1) 

n=0 ^ ' 

which converges for -u G M+, s G C, with either {\z\ < 1, 5J(s) > 0), or {\z\ = 1, ^{s) > 1). 
Special cases include the analytic continuations of the Riemann zeta function and the Hurwitz 
zeta function, (valid for 3?(s) > 1), 

the alternating zeta function (also known as Dirichlet's eta function), and the the Dirichlet beta 
function, (valid for Re{s) > 0), 



00 



k-l 



k=l 



k 



the Polylogarithm, and the Lerch zeta function 

00 f. 

Li„(z) = = z$(^,n,l), L(A,q;,s) = $(exp(27riA),s,Q;). (5) 

k=l 

In the following we will be particularly interested in the evaluation of Lerch transcendent 
on the critical line, when s = 1/2 + it, and for this case \z\ < 1. Thus the special cases of (2), 
do not converge while those of (3), do converge. It is known that the Lerch transcendent can be 
expressed as a Mellin transform: if (|2;| < 1, 9?(s) > 0), or (z — 1, 3?(s) > 1), then 

poo ^(u-l)t 

A(z,s,u) = r(s)^z,s,u)^ — t'-^dt, (6) 

Jo e* - 2; 

and making use of the relation C*(s) = (1 — 2^^*)(^(s), we have the very interesting special 
case, where the Riemann Zeta function can be expressed as a Mellin transform, 

(l-2^-)C(.) = r(«) = $(-1,^,1) = ^/ -r^t'^'dt, (7) 
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In summary, we proceed below to uncover a new unitary change of basis in the case of 
quantum mechanics on R"*", which corresponds to the evaluation of the Mellin transform of a 
wavefunction when s = 1/2+it. This transform corresponds to moving from the (,t G IR+), 
basis to the (p,, e M), basis where we call the latter the Hyperbolic momentum basis. We 
then deduce suitable quantum mechanical potentials and boundary conditions that yield ground 
state wavefunctions which, when viewed in the Hyperbolic momentum basis, give the left hand 
side of (6), when s = 1/2 + it which includes the Riemann-Zeta function as a special case. For 
the Riemann-Zeta function case we find that the relevant wavefunction V^c('^)- x G I^^^ is quite 
simple and corresponds to the ground (bound) state of a deceptively simple ID potential. Such a 
potential might be physically engineered in a variety of systems, e.g. neutral atom mirror traps, 
trapped ion systems. 

As a second case we can also find unbounded, wavefunctions in the \x) basis which, when 
viewed in the basis, give the Riemann Zeta function C(l/2 — iPri), precisely. In this case 
both the original and transformed states are unbounded. However as the basis transform is 
unitary, and thus preserves inner products, we would conjecture that the analytic continuation 
of these states gives results identical to the analytic continuation of the Riemann Zeta function 
onto the critical line. 

Finally we show how one can obtain bounded wavefunctions which can approximate the 
Riemann Zeta function in a more systematic (though somewhat complicated) fashion. We fi- 
nally present a scheme based on the Degenerate Parametric Oscillator to measure the Hyperbolic 
momentum in a physical system corresponding to the Degenerate Parametric Oscillator. 

One interesting result of this work is that we are able to find a physical quantum system 
whose ground state wavefunction, when viewed in the appropriate basis, possesses zeros whose 
locations exactly match those of the Riemann Zeta function on the critical line. This tantalising 
result can be compared with the Hilbert-Polya formulation of the Riemann Hypothesis: to find 
a self -adjoint linear Hermitian operator whose spectrum matches the locations of the zeros of 
the Riemann Zeta function on the critical line (4). Evidence for the latter seems quite strong 
following extensive numerical analysis of the spacings between adjacent zeros of the Riemann 
Zeta function and their comparison with the statistics of a Generalized Unitary Ensemble (GUE 
- or eigenvalues of random Hermitian matrices) (4-8). It is not clear to us how our construction 
has relevance to the Hilbert-Polya formulation of the Riemann-Zeta hypothesis but the existence 
of relatively "simple" quantum wavefunctions with randomly distributed zeros seems intriguing. 

2 Quantum Mechanics on 

We now turn to the description of quantum mechanics on the half-line x E There have 
been numerous expositions of the quantum mechanics of a particle moving in M+. Almost 
without exception, however, these analyses are executed in the position basis |a;). Here, we 
instead seek to make a closer correspondence between the conjugate phase-space operators 
X, p, [x,p\ — ih, for a particle moving on the full-line x e R, with analogous conjugate phase- 
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space operators on the half -line. As the latter must represent physical observables they must 
be represented by Hermitian Self- Adjoint operators. To determine whether an operator A, is 
Hermitian one must examine the action of the operator A, and its adjoint A\ but also the 
domains (■'tp(x) G ^^^[0, cxo], if \^(x)\'^dx is finite), of these operators, D{A), and D(A^). 
If the actions of A and are identical, i.e. (AV,^) = {4>,Ailj), V0,t/' e D{A), then the 
operator A, is Hermitian provided that D{A) e D{A'). Further it is Self-Adjoint if the actions 
match and so too do the domains of the operator and its adjoint, -D(A^) = D{A). An operator 
can be Hermitian but not Self- Adjoint and this will give rise to unphysical situations in that the 
action of may change the domain of the wavefunction, i.e. to different boundary conditions 
etc. 

Here we follow closely the work of (9), and (10), to find that the obvious candidate for the 
momentum operator ~ —ihd/dx, conjugate to a; e R"*", is not Hermitian and Self -Adjoint. 
Moreover, there is no self-adjoint extension possible for this operator. This is to be contrasted 
to the case of a particle confined to a finite interval x E [a, b] (11, 12), or on the entire real line 
with delta-function potentials or with the origin removed (13, 14). We instead proceed to define 
a new momentum operator p^, the "Hyperbolic momentum" which corresponds to dilations of 
the X coordinate or displacements of the "Hyperbolic position" rj = In x, coordinate. 

Before examining the case of p^, let us examine the somewhat easier case of the kinetic 
energy operator on the half -line, 

h? d"^ 

° 2m dx'^ ^ ^ 

As mentioned above we have to pay particular attention to the domain of operators on which 

Hq is allowed to act. We define 'iIj{x) E Vt = C2°([0, cxd]), to be the class of functions which 
are infinitely differentiable (absolutely continuous, i.e. f{x) = ^^{df / dx)dx + /(O)), and have 
finite norm on the interval [0, oo\. If we now consider two functions 0i, 02 £ ^> we can 
examine the Hemiticity of Hq, by integrating by parts. 



(//oV2,0l)-(02,Wl) = 0^(0) 



dx 

0;(O)0i(O) 



x=0 

0'2- 



^1 



0i(O) , (9) 
(10) 

x=0 



If one chooses boundary conditions (A) Ca = 02(0) = 0i(O) = 0, or (B) = 0i(O) = 
d(j)\/dx\Q = 0, then Hq is Hermitian. However in case (A) the domains of Hq and H^ are 
fixed to be identical by C^, and thus Hq is also self-adjoint. In case (B) D{Hq) is reduced to 
01 e Q n Cb, while D{hI) = fl, and thus the operator Hq is not self-adjoint. From (10), 
a more general class of boundary conditions which give an Hermitian and self-adjoint Hq, is 
when — K E R. This choice of boundary condition extends the operator Hq, into a 
one-parameter class of operators and this process is known as a self-adjoint extension. 

If we now repeat the above analysis for the operator p^; = —ihd/dx, again assuming 0i, 02 e 
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Q, we have 



{ci);'cj)i + cj);(j)[)dx , (11) 
= M^*2<l>iL-^*2<l>i\o] ■ (12) 

The first term in (12) vanishes due to the nature of Q, and thus for the second term to vanish 
one must set either 0i(O) = 0, or 02(0) = 0. In either case the resulting operator will be 
Hermitian but will not be self-adjoint. It is not obvious that we can choose alternative boundary 
conditions that will yield a self-adjoint extension of this non-self-adjoint operator. Fortunately 
there is a theorem of von Neumann (14, 15), which tells us whether such an extension exists 
or not. The method is quite simple and involves comparing the dimensions of the so-called, 
"deficiency subspaces", A/±, of the operator in question, e.g. A. These spaces are defined by 

Af±^{iPeD{A^), A^iP^z^iP, ^{z±)^0] , (13) 

and to find (n+, n_), the dimensions of these spaces, one just looks for the number of normal- 
isable independent solutions to A^ijj = ±.i'^ip, for 7 real and positive. For p^. = —ihd/dx, we 
easily have t/j± ~ exp(±7a;/ h), and for x e R+, n_) = (0, 1). From (15), if the deficiency 
indices are different then the operator is not self-adjoint and moreover no self-adjoint extension 
is possible. If we instead had considered the case of the whole real line then nJ) = (0, 0), 
i.e. neither wavefunction is normalisable on the entire real line, and in this case not only is Px, 
self-adjoint, it is termed essentially self-adjoint as its domain is the entire G L^[— 00, -|-oo]. 



3 Hyperbolic Coordinates: Dilatons Displacements 

That the typical expression for the momentum operator p^ — —iM/ dx, does not correspond to 
a physical operator on M+, seems to have only generated little attention in the literature. On the 
other hand, as [x, px] — ih, the action of px on x, yields a displacement, 

eiap./h^^-iap./h ^ ^ ^ ^ , (14) 

and this cannot represent a physical operation as it alters the domain x e M"*". Instead we 
examine the operator which generates dilations of the x operator: 



where the last line is evaluated in the x representation. This operator, as a generator of spatial 
dilations has already been studied as part of a larger space-time conformal transformation in 





.h 


= XPx - 


- 1— 
2 


— —ih 


d 




dx 



+ 



(15) 
(16) 
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certain non-relativistic quantum mechanical problems by de Alfaro, Fubini and Furlan {16), and 
Jackiw {17). The operator also appears in work by van Winter {18), and Berry and Keating {19), 
where they take this operator to be the Hamiltonian of a system. Indeed, although connections 
between the Riemann Zeta function and eigenfunctions of the operator H = xpx were made 
in {19, 20), our treatment differs significantly in that we consider the operator = xpx + PxX, 
to be a new conjugate variable describing dilation/conformal momentum. From [x, p^] = ihx, 
the Sack algebra {21), we can see that the action of this new momentum on x, is as expected: 

^ ^g/. . (17) 



To test for Hermiticity we examine 

i^=(pl01,02)-(01,P,02) , (18) 

and by setting (f)i — U/ -^/x, 02 — V/ we find 

K = C/*y|~ , (19) 

which vanishes if V{Q) = f/(0) = 0, and thus is Hermitian for ~ f/ {x) / y^, {(6 G = 
L^[M"'", dx, U (0) = 0]}. To see if the operator is self-adjoint we use the von Neumann test and 
examine p^-^ = ±i7'0. We find -0 ~ V^/ \/x, but V{x) = A±x^^/^, neither of which are in Q, 
and thus we have the deficiency indices nJ) — (0, 0), indicating that p^ is an essentially 
self-adjoint operator. Defining i) = Inx, (this will be made more precise later on), we can 
recover the standard Heisenberg algebra and displacement operation, albeit in the "exponent 
space": [r),p,,] = ih, and 

^iap^/hf^^-iapr,/h ^ ^ ^ _ ^20) 

Thus by making the unitary transformation of the quantum mechanics on the half line de- 
scribed by the "conjugate operators", (£,p,,), to the exponential representation on the full-line, 
(In X, pr,), we can regain the familiar Heisenberg algebra. Dilations of x now become displace- 
ments of ?7 = In X, and one can formulate all of the standard QM on ]R+, in this Hyperbolic 
phase- space. 

We now evaluate the transition matrix elements {Pr,\x), {Pnlv): {Pvl'^)' resolutions of 
unity in the new basis in order to obtain the correct measures. We will also obtain a displacement 
operator £)( A, //), in the (r^, p^), phase-space and from this a Wigner function pseudo-probability 
representation of a wavefunction in this phase-space. Armed with these tools we will then be 
in a position to seek wavefunctions where (pr,|0) are related to the Lerch transcendent (and 
Riemann Zeta function in particular), on the critical line. First we look at the eigenstates of p^ 
by taking (a;|p^|p^) = p,,(a;|p^), and where ^'^''(a;) = VP'^{x)/^/x, to obtain 

{x\pr,) = -^x'P'^/'^-'^ . (21) 
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As [fj, pr,] — ih, we can define the eigenstates of the coordinate conjugate to pr, via the Fourier 
transform 



1 ^+00 

l^) = ^y^ dp.e-'^^^/'^lp,) , (22) 



which gives 



1 /• + OO „ipr,/h 



e 



-00 

-u/2 f+00 



2nh 



/+00 
dpr, e-^f''(^-")/'^ = e-"/2(5(77 - u) ^ ^/^5{x - e'') , (24) 
-00 



where we have set u = \nx, and have used 5{u)/ g'{u) = 5{x), x = g{u). We see that the i) 
eigenstates are not merely rescaUngs of the x eigenstates but also include a square-root weight- 
ing. From (24), we see that 

{ri\ip)^ cix^(x)V^(5(x-e'') = 6^/^(2; = 6") . (25) 
Jo 

Let us define the following ket so as to absorb the square-root factor: 

\^ =Vd^\x = e'^)^ , (26) 



where we have inserted the subscript x to emphasise that the ket in question is evaluated in the 
X basis but at the value where x — exp{r]). Using this notation, (21), and the properties of the 
5— function we can show 

/+00 
dv{pl\\v){v\\p') , (27) 
-00 

= S{pI-pI) . (28) 
Thus I = Wjildrj, while from (24), we have 

{m\m) = (/ dx \x) {x\^ \r}2) , (29) 

/•oo 

= / dxe'''/^+''''^^6{x - e'")6{x - e''^) , (30) 
Jo 

= e''i/2+''2/2(5(e''i -e^2) , (31) 
= S{vi-V2) . (32) 

One can then easily show that the new coordinate kets I77), (where now rj e K), are completely 
analogous to \x), on the entire real line, and the dilatonic momentum p^, (which generates a 
dilation of x), generates displacements of rj: 
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and now we can make the definition of 77, (which we previously had as 77 ~ In x), more precise: 



+CX3 



77= / dr)r]\ri){ri\ . (34) 



00 



Armed with the conjugate exponential operators, f), and p^j, where [t], p^J = ih, and the transition 
matrix elements {r]\x), {pri\x), we can take any quantum state represented in the x coordinates 
and represent it in the dilation coordinates: (a;i|p|a;2) {vi\p\V2)- 

4 Mellin transforms and the Lerch Transcendent 

As we mentioned above, one is typically familiar with the Fourier transform appearing in quan- 
tum mechanics representing the unitary change in basis |a;) — > x,p E R. To the best of 

our knowledge the Mellin transform has never appeared in the literature playing a similar role. 
However, as we will see, the above change in basis \x) yields a Mellin transform. The 

Mellin transform and its inverse, of a function / is given by 

{Mf}{s) = 4>{s)= / x^f{x)- , (35) 

Jo X 

{M-V} {x) = f{x) - — / x-'(t>{s)dx . (36) 

The Mellin transform may also be defined via the Fourier transform via, 

{W}(s) = {^/(e-^)}Hs) . (37) 
From (21), we see that the representation of a given wavefunction {x\ip), in the basis is, 

POO 

{Pr,\lp) = / dx {prj\x) {x\lp) , (38) 
^0 

= ^= / dx7jj{x)x-^/^~'P^ , (39) 
v27r Jo 



^ {M^l^} {s^ 1/2 -ipr,) . (40) 



^27r 

We now return to our previous expression for the Lerch transcendent (6), we see that 

E{z,s = 1/2 -ipr„u) = r{s = 1/2 -ip^)^z,s = 1/2 -ip^,u) (41) 

roc _(u_i)t 

= / dt— r^/^-iPv , (42) 

Jo e* - z 
= {Mx}{z,s^ 1/2 -ip„,u) , (43) 
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where 

g-(n-l)t 

Xiz,t,u)^— . (44) 

Thus, to simulate the three parameter family of Lerch transcendents 1/2 — ip^, u), we take 
the Mellin transform of xi^: t-, u), evaluated on the critical line s = 1/2 — ip^. However, by 
(40), this can be viewed as the representation of a wavefunction t/j, where 

g— Ma; 

{x\ijj(z, u)) = i/j.^uix) = N(z, u)- — , (45) 

L — ze 

with the normalization factor N{z, u), and where x E M"*". We can now ask about the quantum 
system where the wavefunction (45), could appear. More specifically, can we find a potential 
V{z, u; x), for which (45) is an eigenstate of the half-line Schroedinger equation. 



2 dx"^ 



V{x) 



As we have noted in section (2), the Hamiltonian/Schroedinger operator is only Hermitian and 
self-adjoint if one restricts the domain of the wavefunctions to satisfy '0z,u/'0z,uU=o = k e R. 
From (45), we can find 

lim ; . (47) 

Inserting (45), into (46), and letting R — ze"^ jiX — ze~^), we can solve for the potential 
function, 

V{x) ^E + ^[{u + Ry + R{R + 1)] . (48) 

From the form of (48), we can add a constant to E, without changing the condition that (45) 
solves (46). We add —{E + so that limj;_>oo V{x) — 0, and thus we can set 

V{x)^+^[{u + Rf -u'' + R{R+l)] , (49) 

where now E = —u'^/2, and liuij.^^ V{x) = 0. We note that 

N 1 uz 1 z(l + z) 

and this becomes singular when z — 1. This singularity at 2; = 1 is to be expected as our 
analogy breaks down since (6) holds for z — 1, only when 3?(s) > 1, but above we assumed 

s — 1/2 — iprj. 

A particularly interesting case is when z = —1, u = 1, as the Mellin transform of the 
wavefunction 

(^I^c)=A/'y^ , (51) 
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where the normaUsation factor J\f ^ = —1/2 + In 2, gives via (7), the Riemann Zeta function 

(P.l^c) = A/"!! - 2i/2+''^'")r(l/2 - ^p,)C(l/2 - iPr,)/V2i . (52) 

However the amplitude of this wavefunction, now represented in the Hyperbolic momentum 
coordinate p^, falls off rapidly with p^, due to the F-function in (52). From (50), we see that 
{x\il),^),\s an eigenf unction of the potential 



with the boundary condition 



lim 



+ 1 



1 ^ 
tann — 



X e 



1 

2 



(53) 



(54) 



+0+ ip^ix) 

The prefactors appearing in (52), possesses no zeros and thus the location of the zeros of the 
wavefunction (p^jV'c), corresponds exactly with those of the Riemann Zeta function. However 
these prefactors, (1 — 2^^*)F(s), when s = 1/2 — ip^i, damp as exp(— p^). Thus the resulting 
wavefunction is mostly peaked around p,, ~ 0, and the finer details of the zeros of the wavefunc- 
tion are very difficult to observe for large p,,,. For clarity, since (pr?|^) ~ C(l/2 — iPri) e^^'', in 
order to distinguish the zeros of the wavefunction we plot in Fig. 1, the magnitude of the unnor- 
malised wavefunction (in a logarithmic scale), rescaled by p^. We see perfect correspondence 
between the wavefunction and the Zeta function zeros. Below we will show other forms for 
ip({x), which yield much closer approximations to the Zeta function, i.e. where the prefactors 
vary more slowly with the Hyperbolic momentum p^. 



5 Exact quantum simulation of the ( function 

From Fig. 1, we see good correspondence between ^^(p^), and the Riemann Zeta function on 
the critical line. However the correspondence is detracted by the exponential damping given 
by the F-function appearing in (7). It would be preferable if one could find a wavefunction in 
the a; —representation {x^tlj), whose representation in the Hyperbolic momentum (p,,|^), gives 
a closer match to the Riemann Zeta function than afforded by (7). Surprisingly a wavefunction 
which yields a perfect match can be found. To describe this wavefunction we make reference to 
a set of qubit basis states proposed by Gottesman, Kitaev and Preskill (GKP), (22), for quantum 
information processing using continuous variables. On the entire line, x e IR these states are 
defined by 



+ 00 



|0) = J2 = = J2 ne-^n) , (55) 

n=— oo n=— oo 
+00 oo 

py = ^ |x = 2en + e)= ^ = Trr^n) , (56) 
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9 e [0, 1], which are infinite superpositions of infinitely squeezed states (i.e. eigenstates of x, 
andp), which are unnormalised and also possess infinite energy. Such continuous variable states 

have been proposed to encode a qubit \q) = a|0) + 6| 1) (22), and recently an error analysis {23), 
and various physical methods for their synthesis (24-29), have appeared in the literature. As 
we are concerned with quantum mechanics on the half-line we set 6' = 1/2, and only consider 

the 5-function sum in |0), over a strictly positive domain, 

+ 00 

|0) ^J2\x = s) . (57) 

n=l 

Following from (21), we can easily find 

n=l 

where s — 1/2 — ip^^, which is precisely the Riemann Zeta function evaluated on the critical 

line. However, just as the state |0), in the x— representation is unbounded and possesses infinite 
energy, the Riemann Zeta function for 3?(s) < 1 is also formally unbounded and is defined 
there via analytic continuation. Thus although formally we have a wavefunction whose p^- 
representation is precisely the Zeta function, the unboundedness poses severe difficulties in its 
physical interpretation. 



6 Intermediate simulation 

We have seen a wavefunction whose -transform yields the Riemann Zeta function up to a 
rapidly damped multiplicative function and another wavefunction whose p^-transform yields the 
Riemann-Zeta function precisely. We now instead search for an intermediate wavefunction in 
the representation which does not possess any pathologies such as unbounded norm, infinite 
energy etc and which gives the Zeta function up to some multiplicative factor which damps far 
slower than above, i.e. slower than r(l/2 — ip^) ~ exp(— p^). As we shall see, there are an 
infinite number of such examples. For the most part however, we will be unable to give simple 
analytic formulae for these wavef unctions. 

To find these wavefunctions we return to Eqn. (7), 

/•oo 

(1 - 2^-)C(.)r(.) = -r^t'-'dt . (59) 

This is a special case of a more general Mellin transform formula (30), 

^ IE (-l)""V(^ri) ] 1 (s) = (1 - 2^-) {Mf} {s)C{s) , (60) 



,n=l 
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where f{x) is any function. We see that the particular case of (59), is obtained when we choose 



fix) = exp(-x) , {Me--} (s) = r(s) . (61) 

We thus seek a function (say), g, whose Mellin transform {A4s}{s), does not decay rapidly 
with t, for s — 1/2 + it. An example of such a function is the one parameter family of functions 

,x 1 +xcos((h) 
s(^-^)=m^) , +2.00^1.^ ■ (62) 

where N{(j)) is a normalisation factor. This function has the Mellin transform (30), 

E{s,<f>) = {Mgm (63) 

N{(f))7r 



2 sin(7rs) 



[(cos0 - i| sin0|)^ + (cos0 + i| sin0|)*] . (64) 



The sum. 



E(a;,0)^-5^(-l)V(na;) , (65) 

n=l 

cannot be easily expressed in analytic form and must represent a normalisable wavefunction on 
X G M"*", and thus we choose the normalisation factor N{(j)), such that \T.\'^dx = 1. We plot 
the behaviour of (the unnormalised), T.{x, 0), in Fig. 2. As ^ vr, T.{x, 0), develops an infinite 
number of oscillations within the domain x G [0, 1], but in the following we will only consider 
G [0, 3], and here the function g{x, 0), and its Mellin transform (p), are well behaved. We 
see that the Mellin transform damps quickly with increasing t, for low values of G [0, 2.5], 
but as ^ 3, this function damps more slowly. We expect that by increasing the value of the 
parameter tt, we can extend the Mellin transform's E{l/2 + it, cf)), extent to arbitrary large 
t-values at the expense of using highly oscillatory wave-functions S(a;, 0). If we now follow this 
argument and choose a wavefunction in the x— domain to be {x\x) = N{(f) = 3)S(a;, (p = 3), 
we can move to the p^— representation to find 

(p,|x) = iV(3)(l-2V2+^^"')C(l/2-^p,)S(l/2-zp„3) , (66) 

and we plot vs. — p^, in Fig. 1 (11). Compared with Fig. 1 (I), we can clearly see 

that the wavefunction now has greater extent over the —p^ coordinate and more details of the 
zeros/nodes of this wavefunction are apparent. 

Thus although we have shown that one can find legitimate normalisable wavefunctions on 
the X G R+, axis whose -representation gives closer and closer approximations to the Rie- 
mann Zeta function, these wavefunctions become more and more pathologic e.g. ~ sin(l/x). 



12 



7 Wigner Functions in the (77, p^^) phase space 



The Wigner quasi-probability distribution function W{ri,pri), is a standard tool for vizualising 
the quantum aspects of a quantum state, i.e. when W assumes negative values. The W function 
however, is typically defined on the phase space {x, p), where x, p e R. For quantum mechanics 
on the half-line x e we have shown above that one cannot define a self-adjoint momentum 
operator and thus we cannot construct an associated Wigner function. However in the Hyper- 
bolic representation {r],Pri), where now both e M, and [?7,p^] = ih, we can (as in (9)), 
define a Wigner function to be 

WD{r],p^) = — dr]'{r] + -r]'\p\r] - -7]') e^"'*'" . (67) 

For the two quantum states (52), and (66), we have plotted the associated Hyperbolic phase- 
space Wigner function in Figs. 3-6. Although zeros of the wavefunction do not normally 
correspond to zeros of the Wigner function, from Figure 4 and Figure 6, some correspondence 
seems present. Moreover the p^-marginal probability, |V'^(Pr/)P = / W^(^,P?j)<^^> exhibits the 
engineered Riemann-Zeta zeros. 



8 Physical Realisation 

Above we have given a mathematical derivation of a new type of unitary transformation which 
transforms quantum mechanics on the positive half line to a Hyperbolic representation and in 
the process, executes a Mellin transform '\\){x) {M.'i\}\ (s = 1/2 — ip^^l\ph\. By choosing 
'4}{x), carefully (such that it is the eigenstate of a particular potential V{x), (53), subject to the 
boundary condition (54)), the transformed wavefunction ■ilj{pri), has a nodal structure identical 
to the Riemann Zeta function on the critical line. The boundary condition (54), is not typi- 
cal. However as shown in (31, 32), one can closely approximate such a boundary condition 
by a short-range potential Vbndryix), near the origin. For the boundary condition (54) to hold 
independent of energy, Vbndry, limits to a hard-wall infinite potential with infinitesimal struc- 
ture (52) which enforces (54). Possible physical systems which might admit the engineering 
of such short-range boundary potentials are multi-frequency evanescent-wave mirrors used for 
trapping and reflecting ultra-cold neutral atoms (33). A more challenging task is the physi- 
cal realisation of the measurement of Hyperbolic phase (pr;|V^)- We propose using an indirect 
measurement scheme (34, 35), where the system we wish to measure is coupled to a quan- 
tum probe via an interaction Hamiltonian Hj such that the Hyperbolic phase of the system 
drives a linear displacement of the probe's quantum state. Since [rj,Pr,] — 1, we can formu- 
late Hyperbolic creation and annihilation operators r] = {d + d))l% Prj = {d — a^)/2i, with 
[a, a^] = 1. Taking the probe as a bosonic field described by operators b, b\ [b.b^ = 1 
(i.e. on the full line rather than the half line and these quadrature operators Xpb, Ppb, are now 
well defined self-adjoint operators), then the indirect measurement generated by the interaction 
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Hf — x(a6^+a^6) = 2ix(xp6Pf,—pp6r7), will drive displacements ofthe probe field generated by 
the system's Hyperbolic quadratures. Unfortunately, although (from (16)), = [xp^ +pxx\/'2, 
the Hyperbolic position is effectively the logarithm of the position of the system on the half line, 
•q ~ Inf, and since Inx (— — a)Y(ia*), a; < 2a, the interaction Hamiltonian 

Hf involves coupling the probe to a highly nonlinear function (involving all powers of x), of 
the system's position. Instead we consider the alternative coupling Hamiltonian 



B 



2X [Xph{.xl -pl) +pph{,Xsps+PsXs)\ 
2X [Xph{xl - Pi) + 2ppbPrj] 



(68) 



where x. 



a, a 



\ now describes a bosonic mode defined on the entire real line. We note 



however that (68) acts disjointly on the two half-line "super-selection" sectors x G IR+ and 
a; e M~. From this we see that the probe's momentum is coupled to the Hyperbolic momentum 
in each super-selection sector but the probe's position is no longer coupled to the Hyperbolic 
position. Moreover, (68), is the Hamiltonian for the Degenerate Parametric Oscillator, and 
has been studied by many (36^0), though mostly in the case where the probe field is treated 
semiclassically. To get a feeling for how (68), provides us with an indirect measurement of the 
Hyperbolic momentum of the system mode we examine the quantum Heisenberg equations of 
motion and their semiclassical c-number approximations. 



da 



da 
'dt 



-2ixPa 



*2 



db . ^„ d(3 .2 



(69) 
(70) 



given the initial probe state to be the vacuum \ijjpi,) = |0), and the system state only within 
]R+, i.e. {x\ipsys) = Vx < 0. Setting r = Axt, quadratures for the probe and system 

2(3 = Xpb + ippb, 2a — Xg + ips, and introducing the new variables 2u = x1 — pi, 2w = 
+ f's! ^ = ^sPs, where now v is the semiclassical analogue of the Hyperbolic mometum of 
the system, equations (69,70) can be recast as 



d_ 
d^ 







' 


1/2 








" 




•^pb 


V 










-J^pb 
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W 












Ppb 
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u 










Ppb 










u 















-1/2 







_ Ppb 



(71) 



From this one can prove that (ajpj, + p^j,) + w = 2\l3{t)\'^ + \a{t)\^ = is a constant of the 
motion, while in addition Xpb/xpb = ppb/ Ppb- Using this constant of the motion one can arrive 
at coupled nonlinear dynamics for the probe mode alone. 



2x 



pb 



-Xpb[K - {xlb + pjb)] , 2ypb = -Ppb[K - {xlb + pjf,)] 



(72) 
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with the initial conditions Xpb — Ppb — 0, 2xpb — XsPs, '^Ppb — ^1 — pi. This is a conser- 
vative system of nonlinear ordinary differential equations which corresponds precisely to the 
kinematics of a particle moving in two spatial dimensions {xpb, Ppb), under the central potential 



V{xpb,Ppb) = - ^r^ 



(73) 



pb Ppb- dynamics can be exactly solved via Elliptic functions. One can find 



where = x 
a series solution for small times to be 

1 



Xpb{r) 



-XsPsT 



Ppb{r) = --{x1-pI)t 



(74) 
(75) 



Qualitatively the particle moves away from the origin initially with constant velocity (as near 
the origin the force vanishes), but after a time t* ^ 



4wit=o) Ppb asymptote to 

fixed values as the particle comes to rest at the maximum of the potential V. The initial energy 
of the system w{t — 0), drives the probe up the potential and in the long time limit all of the 
system's energy is expended into displacing the probe. However the dependence of the probe's 
displacement on the initial system parameters, u and v is only linear in the small time limit 
t -C t*. To summarise, from the semiclassical Heisenberg equations of motion we have good 
reason to believe that for short interaction times the system's Hyperbolic momentum via, Hf, 
induces a linear response of the Xpb quadrature. Although Degenerate Parametric Oscillators 
are studied in optical settings it is also possible to engineer (68), between the vibrational modes 
of an ion in an ion trap (38). 



9 Conclusion 

We have found a new type of unitary transformation for quantum mechanics on the half -line 
which transforms one into the Hyperbolic representation. The effect of this transformation is to 
execute a Mellin transform on the wavefunction. Although this work has no immediate impact 
on the Riemann-Zeta hypothesis, it is interesting that we are able to present bounded quantum 
wavefunctions whose nodal properties match exactly those of the Riemann Zeta function on the 
critical line. It is also possible that just as the quantum Fourier transform has found numerous 
applications in quantum computation, this new quantum Mellin transform might lead to new 
quantum applications and algorithms. 
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Figure 1: Comparison of \({s = 1/2 — iprj)\ (blue), and |V'c(Pr?)l (red) wavefunction in (I) 
from Eq.(52) and (II) from Eq.(66). In (I) we plot [log |V'c(Pr?)l]/(~Pr?) (red), to overcome the 
exponential supression of (52), generated by the r(l/2 — ip^) function while in (IIA) we plot 
li'ciPc)] (red) on a linear vs. linear scale and in (IIB) we repeat on a log vs. linear scale. In 
(IIA) and (IIB), the details of the zeros of the wavefunction are more clearly visible for a larger 
range of p^. 
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Figure 2: Graphs of (A) the sum 0) = Yl'^=i fif^^)^ from Eq. (65), and (B) the 
Mellin transform of /(x), 0)| from Eq. (64), for s = 1/2 + it. For low values of ~ 2, 
0) possesses very little structure while for ~ 3, this function develops more complex 
oscillations. For low values of ~ 0, the Mellin transform H is highly damped function of t, 
while for ~ 3, the transform damps much more slowly, thus allowing extended explorations 
of the Zeta function on the — axis. 
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Figure 3: Hyperbolic phase space Wigner function W{r],p,j), of the Riemannn-Zeta wavefunc- 
tion in Eq. (52). 
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(A) (B) (C) (D) 




Figure 4: Details of the Hyperbolic phase space Wigner function W{ri, p^i) of the state (52). (A) 
W{r],pr)), (B) In \W{r],pr))\, (C) Pr, vs. / W{r],pr,)dr], and (D) \({l/2 - ipr,)\. One can see faint 
traces of the Riemannn-Zeta zeros of the wavefunction within the Wigner function. 
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r] 



Figure 5: Hyperbolic phase space Wigner function W{r],p,j), of the Riemannn-Zeta wavefunc- 
tion in Eq.(66). 
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Figure 6: Details of the Hyperbolic phase space Wigner function W{i], pr,) of the state (66). (A) 
W{r],p^), (B) In \W{r], p. (C) p^ vs. f W{r],p^)dr], and (D) |C(l/2 - ipr))\- As compared 
with Fig. 4, the Riemannn-Zeta zeros are more clearly displayed in this wavefunction. 
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